Field-induced Bose-Einstein Condensation of triplons up to 8 K in Sr 3 Cr 2 08 
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Single crystals of the spin dimer system SraC^Os have been grown for the first time. Magne- 
tization, heat capacity, and magnetocaloric effect data up to 65 T reveal magnetic order between 
applied fields of H e i ~ 30.4 T and H c2 ~ 62 T. This field-induced order persists up to T™ ax ~8Kat 
H ~ 44 T, the highest observed in any quantum magnet where H C 2 is experimentally-accessible. We 
fit the temperature-field phase diagram boundary close to H c i using the expression T c = A(H-H c i)". 
The exponent v = 0.65(2), obtained at temperatures much smaller than T™ ax , is that of the 3D 
Bose-Einstein condensate (BEC) universality class. This finding strongly suggests that Sr3Cr20s is 
a new realization of a triplon BEC where the universal regimes corresponding to both H c i and H C 2 
are accessible at 4 He temperatures. 

PACS numbers: 73.43.Nq, 75.30.Kz, 75.30.Sg, 75.40. Cx 



The Bose-Einstein condensate (BEC) is a unique state 
of matter in which a macroscopic number of bosons col- 
lapse into a global (identical) quantum statejj, Q. Al- 
though supcrfluid 4 He [3] is possibly the most studied 
BEC system, ultracold Rb gas proved to be the first true 
experimental realization of this exotic state of matter [![. 
Since the original discovery was made, other BECs have 
been identified in optical lattices, superconductors 0], 
and optically-pumped magnon gases @ where the phe- 
nomenon of BEC can be studied at room temperature. 

Giamarchi and Tsvelik first showed that the Hamilto- 
nians of quantum antiferromagnets and BECs were di- 
rectly related by a mapping transformation 6] , follow- 
ing up on earlier work by Matsubara and MatsudaQ- 
Although real lattices always have imperfections and 
symmetry-breaking terms, in some materials these can 
be neglected since they only affect their properties below 
the experimentally-accessible window. This then sug- 
gests that a very good approximation to a BEC can be 
achieved in solid-state magnetic systems. In particular, 
spin dimer systems were perceived to be good candidates 
to display BEC behaviour. These systems exhibit spin- 
singlet ground states with a spin gap to the first excited 
triplet state. This gap can be closed by applying a large 
enough magnetic field, resulting in the generation of a 
macroscopic number of triplet excitations (triplons). If 
the kinetic energy terms dominate the potential energy 
terms in the U(l) invariant Hamiltonian of such a system, 
the resulting ground state can be described as a BEC. 

TICUCI3 was the first spin dimer system to show evi- 
dence for a triplon BEC ground statejl, and this has 



stimulated a flurry of research in triplon BECs[10|, [II 
12 . 1 13 . Il4|. Currently, the primary goal is to search 



for systems displaying BEC behaviour in the space of 
externally-controlled parameters completely accessible in 
the laboratory. This allows the associated quantum 
phase transitions to be studied experimentally. 

The field-induced order in a triplon BEC is nothing but 
the well-known XY-AFM (easy-plane antiferromagnetic) 
state. Such a system can either be visualized as a lattice 
gradually filling with hard-core triplons, or as an XY- 
AFM where the in-plane component of the canted spins 
varies from zero in the paramagnetic state, to a max- 
imum at approximately midway between H c i and H C 2, 
to zero again in the fully-polarized state. The external 
control parameter in both cases is the applied magnetic 
field, and the order parameter is the number of triplons 
or, equivalently, the magnitude of the staggered moment 
M xy . The (H,T) phase boundary in ideal 3D systems fol- 
lows the expression T c = A(H-H C ) 2 / 3 close to the critical 
fields H c i and H C 2- Symmetry breaking terms can affect 
real systems by changing the universality class from XY- 
AFM (easy plane) to Ising (easy axis), in which case the 
critical exponent changes from 2/3 to 1/2, or by changing 
the dimensionality of the system fljj|. 

In this letter, we present the first evidence for a BEC 
ground state in the spin dimer system S^Cr^Os through 
a combination of high field magnetization (M(H)), heat 
capacity, and magnetocaloric effect (MCE) measure- 
ments. This is one of a series of isostructural spi n dimer 
compounds with the general formula A3M20s[3, E3, El 
where A = Ba or Sr and M = Cr or Mn. These com- 
pounds crystallize in the R-3m space group at room tem- 
perature, and the magnetic M 5+ ions are tetrahedrally- 
coordinated with O 2- ions. The M 5+ ions may carry 
spins of either S = 1/2 or 1, and as a result of having 
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only one nearest neighbour (NN) magnetic ion they form 
spin dimers along the c-axis. While the A 2+ ions are es- 
sentially isolated from the MO 4- tetrahedra, they seem 
to play an important role in controlling the NN M-M dis- 
tance and the value of the associated spin gap. For exam- 
ple, while Ba 3 Cr 2 C>8 has a NN Cr-Cr distance of 3.9765 
A and a spin gap of 15.6(3) Kpjj, Sr 3 Cr 2 08 has a smaller 
NN Cr-Cr distance of 3.842 A and hence a larger spin 
gap of 61.9(1) k[T3|. This is consistent with the expec- 
tation that a smaller NN distance between the magnetic 
ions should yield a larger antiferromagnetic intradimer 
exchange interaction. The stronger intradimer exchange 
interaction in the Sr-Cr system pushes the temperature 
scale of the magnetic order to the highest values observed 
yet in a system where H C 2 is experimentally-accessible. 

Single crystals of S^Cr^Os were grown by the optical 
floating zone method at McMaster University. We first 
mixed together stoichiometric quantities of SrC03 and 
Cr203. The resulting powder was shaped into 10 - 12 
cm long rods, which were annealed for 72 hours in air at 
1200°C and then quenched to room temperature. The 
rods cracked during the annealing process, so they were 
remade following this and subsequently annealed for ^5 
hours in high purity Ar gas at 1200°C to strengthen them 
for the floating zone growths. Since S^C^Os is congru- 
ently melting in Ar gas, the crystals were grown at a very 
high rate of 20 mm/h with an Ar pressure of ^170 kPa 
and a temperature of 1470°C. The phase purity of the 
resulting crystals was confirmed by powder x-ray diffrac- 
tion, and they were aligned by Laue x-ray diffraction. 

The anisotropy of the g-factor was measured via elec- 
tron paramagnetic resonance (EPR) at 70 GHz using a 
resonant cavity mounted in a cryogenic goniometer, so as 
to rotate the sample and resonant environment in situ. 
The temperature was maintained at 30 K to ensure suffi- 
cient thermal population of the triplet state using a 4 He 
flow cryostat. The magnetic field was applied using a 
superconducting solenoid and calibrated using the con- 
ventional EPR-marker DPPH. The g-factors were found 
to be g c = 1.938(6) and g ab =1.950(1). 

Fig. [1] depicts dc susceptibility measurements of our 
crystals. At high temperatures, the data follows a Curie- 
Weiss law. Below ~ 40 K, the susceptibility decreases 
sharply with decreasing temperature, which is character- 
istic of non-magnetic spin-singlet ground state systems. 
The susceptibility measurements also reveal that this sys- 
tem is nearly isotropic. We fit the susceptibility data to 
an interacting dimer model of the form: 
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where N^ is Avogadro's number, /is is the Bohr mag- 
neton, Jo is the intradimer exchange constant, and J' 
is the sum of the interdimer exchange constants. The 
last two terms represent contributions from Van Vleck 
paramagnetism/core diamagnetism and impurity/defect 
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FIG. 1: (color online) dc susceptibility measurements of 
Sr3Cr2 08 with an applied field of 1 T. The data was fit to 
an isolated dimer model. The inset shows the magnetization 
M(H) as measured in a 65 T short pulse magnet. 



spins respectively. We fixed the g-factors in the fits to 
the values obtained by EPR. From the fitting, we found 
only ^0.8% free S = 1/2 moments, and estimated the ex- 
change constants to be J = 61.5(3) K and J' =12(2) K. 
Note that the latter are actually best determined by in- 
elastic neutron scattering [20j. 

Magnetization data was collected using an extraction 
magnetometer in conjunction with a 3 He fridge and a 65 
T short pulse (35 ms) magnet at the National High Mag- 
netic Field Laboratory (NHMFL) in Los Alamos, New 
Mexico. The Fig. 1 inset shows M(H) for the H±c ori- 
entation while sweeping the field up. M(H) has a very 
small value up to H c i = 30.4 T due to defects/impurities 
in the sample. Above H c i = 30.4 T, the spin gap is closed 
resulting in a magnetically-ordered regime between H cl 
and H c2 = 62 T where M(H) increases roughly linearly 
with H. Finally, above H C 2 M(H) saturates as the system 
reaches the fully-polarized state. 

To further characterize the field-induced order, heat 
capacity and MCE measurements were performed using 



a home-built calorimeter [12J, [2l|. These measurements 
were conducted in a 3 He fridge placed in the core of a 35 T 
resistive magnet at the NHMFL in Tallahassee, Florida. 
The MCE data was collected for the H_Lc orientation, 
sweeping the field both up and down at 3 T/min. The 
quasi-adiabatic MCE yields a sharp change in the sam- 
ple temperature when entering or leaving a magnetically- 
ordered state to ensure entropy conservation[22|. 

Fig. 2(a) shows a color contour plot of the entropy as 
a function of temperature and field. AS was calculated 
from the MCE data using the equation: 



AS = 



k(T — Tbath) 
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FIG. 2: (color online) (a) Color contour plot of the entropy 
S/T of Sr 3 Cr 2 08, as calculated from the MCE data. The 
phase boundary determined by MCE is shown with black 
dots, (b) Selected heat capacity data as obtained by the dual 
slope method. For 33 T, data was also collected by the stan- 
dard relaxation method (black triangles). 



where k is the thermal conductivity of the thermal link in 
the calorimeter, T is the sample temperature, and T^ at h 
is the bath temperature. We measured function 
of T at both and 33 T and found very little difference 
in the two curves, so the 33 T thermal conductivity data 
was used in all our calculations. The absolute entropy 
was then obtained by a calibration process using 31.5 T 
heat capacity data. We chose to plot S/T to improve 
the contrast since the entropy associated with the transi- 
tion drops quickly with decreasing temperature. As one 
crosses the phase boundary at a constant temperature, a 
peak in the entropy is clearly evident. This results be- 
cause both the low-field gapped phase and the ordered 
phase above the transition should have small entropies. 
However, in the proximity of the transition we likely have 
fluctuations of the order parameter characteristic of a sec- 
ond order transition and hence a much higher entropy. 

Heat capacity measurements are presented in Fig. 2(b) 
for the H_Lc orientation. The majority of the data was 
collected using the dual slope method with a large delta 
t(23|. The main advantage of this technique is that data 
can be collected quickly which is important for efficient 
use of power in high field resistive magnets. However, this 
often comes at a cost of reduced accuracy, and so data 
is generally collected for at least one applied field using 
the standard relaxation technique to compare the two 
methods directly. In the present case, this was done at 



33 T (black triangles in Fig. 2(b)). The excellent agree- 
ment between the techniques suggests that the dual slope 
method estimates the heat capacity well in Sr 3 Cr 2 08. 

The emergence of a lambda-like anomaly was clearly 
observed with an increasing applied field beyond 
H c i ~ 30.4 T; this provides unambiguous evidence for the 
field-induced order in this system. The lambda anomaly 
is small for fields close to H c i, but becomes much more 
prominent with increasing applied field as seen in Fig. 
2(b). Heat capacity data was also collected for the H||c 
orientation. Similar lambda anomalies were observed in 
that case, and the transition temperatures were system- 
atically shifted down by only ~ 0.2 K for a given applied 
field. This indicates that S^Cr^Os is a very isotropic, as 
previously suggested by the susceptibility measurements. 

The spin gap in S^Cr^Os is quite large, and so DC 
fields provided by resistive magnets can access only a 
very small region of the phase diagram. For this reason, 
we also performed MCE measurements in both the 50 T 
mid-pulse (300 ms) and 60 T long-pulse (2 s) magnets in 
Los Alamos. For these measurements, a tiny S^Cr^Os 
single crystal (200 x 100 x 50 um 3 ) was inserted into 
a small, home-built calorimeter [24|. To the best of our 
knowledge, these were the first successful MCE experi- 
ments performed in a pulsed field facility. 

Fig. 3 presents a phase diagram for the H_L c orien- 
tation of S^C^Os as obtained from our specific heat, 
MCE and M(H) experiments. The asymmetric nature 
of the dome is quite unusual for a triplon BEC system; 
further studies are required to understand this feature. 
From the heat capacity measurements, the ordering tem- 
perature for a given field was taken as the peak of the 
lambda-like feature. From the MCE data, the ordering 
temperature for a given field was determined by locating 
the extrema in the first derivative of the sample tempera- 
ture with respect to field) 12. 1 211]. Some selected upsweep 
MCE curves are superimposed on the phase diagram. For 
clarity reasons we show only one downsweep MCE curve; 
the sample temperature changes in the opposite manner 
to the upsweep data near the transition as expected. No 
irreversible/dissipative mechanisms typical for a first or- 
der phase transition were observed. Many MCE scans 
were completed at low temperatures so it would be pos- 
sible to determine a critical exponent. The fitted data 
all must lie in the universal regime for the obtained crit- 
ical exponent to be meaningful, and an earlier Monte 
Carlo simulation study [25| suggested that the universal 
regime extends no higher than ^0AT" lax , where T™ ax is 
the maximum temperature of the magnetically-ordered 
regime. In the present case, we have determined that the 
top of the dome lies around 8 K. To ensure that we stayed 
within the universal regime, only the MCE data below 
2.7 K was fit to a power law of the form: T c = A(H-H c i) 1/ . 

Although one can perform a three parameter fit, H c i 
and v are not independent variables. For this reason, we 
resorted to a series of two parameter fits where H c i was 
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FIG. 3: (color online) S^CnOs phase diagram with selected 
MCE traces superimposed and the ordered region shaded as 
a guide to the eye. Due to small temperature changes at the 
transitions, the MCE traces have been magnified. The inset 
shows the result of fitting the data near H c i to a power law. 



systematically fixed at different values and v was deter- 
mined from the fitting. Our best fit yielded v =0.65(2) 
with a corresponding H c i = 30.40(1) T. The former is in 
agreement with the expected value of 2/3 for a 3D BEC 
universality class, strongly suggesting that S^C^Os is a 
realization of a new triplon BEC. 

Finally, let us comment here on the nature of the 
magnetically-ordered state in S^C^Os- The spin struc- 
tures in the ordered states of the isostructural systems 
BasC^OsdH and Ba 3 Mn 2 08 21 1 have both been deter- 
mined previously, and the results were quite different. 
In the case of BasM^Os, a combination of single ion 
anisotropy and geometric frustration lead to both an in- 
commensurate spiral spin phase and an Ising-type mod- 
ulated structure in the ordered regime 21| . The shape 
and structure of the phase diagram for that material 
depends on the orientation of the system relative to 
the applied field. In sharp contrast, neutron measure- 
ments revealed a commensurate, collinear spin structure 
in BaaC^Osjm. This is likely a result of the geometric 
frustration being relieved due to a structural distortion 
courtesy of a co-operative Jahn- Teller effect |26j]. A simi- 



lar Jahn- Teller effect has been reported for Sr^Cr9Qs[27l|. 
suggesting that the spin structure of the ordered state is 
very similar to the case of BaaC^Os- 

In summary, we have determined the full (H,T) phase 
diagram for S^C^Os from both thermodynamic mea- 
surements and high-field magnetization. This phase di- 
agram provides direct evidence for field-induced order 
in this system between H c i ~ 30.4 T and H C 2 ~ 62 T, 
with an unprecedented maximum ordering temperature 
T™" 1 ~ 8 K. The phase boundary near H c i was fit to 
a power law. Our obtained critical exponent of 0.65(2) 



agrees with the 3D BEC universality class, and so this 
strongly suggests that S^Cr^Os is a realization of a new 
triplon BEC system. Since a large portion of the uni- 
versal regime in this material can be accessed via 3 He 
systems, both the H c i and H C 2 quantum critical regimes 
can be studied in detail. This may help us to understand 
the asymmetric nature of the phase diagram and other 
more general aspects of triplon BEC phenomena. Note 
that no experimental criticality studies near H C 2 have 
been reported for triplon BEC systems to-date. 
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